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CHANGING GEARS: 

ISOSPECTRALITY VIA EIGENDERIVATIVE TRANSPLANTATION 

PETER DOYLE AND PETER HERBRICH 


Abstract. We introduce a new method for constructing isospectral quantum graphs that 
is based on transplanting derivatives of eigenfunctions. We also present simple digraphs with 
the same reversing zeta function, which generalizes the Bartholdi zeta function to digraphs. 


1. Introduction 

Quantum graphs are singular one-dimensional manifolds equipped with self-adjoint oper¬ 
ators, whose spectra allow for explicit computations. Following |BE09] . let G = {V,E,l) 
be a hnite metric graph with vertices V = {ui, U 2 ,..., u„}, edges E = {ci, € 2 , ■ ■ ■ ,em}, and 
edge lengths I = (/i, h, ■ ■ ■, Im) £ • In particular, e* is parameterized by Xj G [0, /j], which 

determines an orientation of the edges. However, the differential operator of interest acts as 
Aj = —^ on (^^(0, li) and is thus invariant under transformations of the form x[ = U — Xj, 
which allows to regard edges as undirected. The self-adjoint extensions of the correspond¬ 
ing symmetric operator with initial domain ©™^RQ(0,/j) can be parameterized in terms of 
boundary conditions at the vertices |KS99a] . which turn the metric graph into a quantum 
graph. Unless otherwise stated, we consider Kirchhoff-Neumann conditions, which require 
that for each vertex Uj, functions on adjacent edges take the same value at Vi, and the sum of 
their outgoing derivatives at Vi vanishes. In particular, leaf vertices carry Neumann bound¬ 
ary conditions. It is well-known that hnite quantum graphs have discrete spectrum |BK13] . 
Quantum graphs are called isospectral if their spectra coincide, including multiplicities. 

Recently, it was discovered |OB12] that the graphs in Figure fTAl have the same set of eigen¬ 
values when viewed as either edge-weighted combinatorial graphs or quantum graphs with 
edge weights or lengths (a, a, 6, 6, c, c), respectively. These graphs hrst appeared in [MM03j 
where a certain line graph construction was applied to the 7i pair of Dirichlet isospectral 
planar domains in |BCDS94] . The second author |Herl5] has generalized this construction to 
manifolds with mixed Dirichlet-Neumann boundary conditions and revealed its connection to 
the graph-theoretic characterization of the famous Sunada method |Sun85] given in |Herll) . 

It is worth mentioning that the main arguments in |OB12) are based on the widespread 
misconception that eigenfunctions on quantum graphs with Kirchhoff-Neumann conditions 
are determined by their values at the vertices |KS99bl ISS061IBS07] . For example, if a, h, 
and c are integer multiples of some r > 0, then each of the graphs in Figures fTAl and [TBl has 
countably many eigenfunctions which are supported on its central cycle of length a + b+c and 
vanish at all vertices. Instead, |OB12] asserts that if /*: [0,/j] —)■ M denotes the restriction 


2000 Mathematics Subject Classification. 34B45, 05C50, 58J53. 

Key words and phrases. Isospectrality, Quantum graphs, Markov chains. Digraphs, Zeta functions. 

We thank Ram Band who conjectured that the graphs in Figure [lAl feature a new kind of transplantation. 































2 


PETER DOYLE AND PETER HERBRICH 







(C) Pi = p'^ + wt'^ and U = p- — t', or vice versa 


Figure 1. Eigenderivative transplantation. At the leaf vertices, we impose 
Nenmann conditions. At the polygon vertices, we reqnire continnity as well as 
Pjih) + = p'fc(O) + wt'fc(O) and +wi'-{lj) = p'^{0) + wtl(O). 


of an eigenfnnction to the edge e* with eigenvalne A > 0, then 

(1) fi{xi) = — / (fi{0) sin(\/A(h - Xi)) + fiih) sin(\/Axi)) . 

sin(vAti) ^ ' 

Similarly, |KS99b| snggests to ignore edges Cj satisfying y/Xk = kn for some fc G when 
determining whether A is an eigenvalne, while snch edges contribnte the bonndary conditions 

(2) Mh) = (-l)Vi(O) and f'ih) = (-l)V'(O). 

In Section [21 we introdnce a method for constrncting isospectral qnantnm graphs that 
avoids explicit compntations and thns bypasses the above-mentioned shortcomings. It is 
a derivative of Bnser’s transplantation method |Bns86] . which itself can be viewed as the 
combinatorial incarnation of the Snnada method jSunSSj. In contrast to the latter, onr 
method can prodnce pairs of qnantnm graphs withont common covers. In Section [3l we 
relate onr method to random walks on combinatorial graphs. In addition, we present related 
pairs of non-regnlar simple digraphs that have the same reversing zeta fnnction as introdnced 
in |Herl4] . Thns, these digraphs exhibit a noteworthy degree of spectral indistingnishability. 

2. Gear graphs and eigenderivative transplantation 

Definition. An n-gear is a qnantnm graph with 2n vertices that is comprised of a polygon 
with n sides of lengths ■ ■ ■ Jn) G as well as n leaf edges of lengths (Zi, / 2 , • • •, G), 

called teeth, snch that the ith tooth is adjacent to the zth polygon side. The dnal n-gear is 
obtained by attaching each tooth at the other vertex of its corresponding polygon side. 

Fignres [TA] and ITBI show pairs of mntnally dnal 3-gears with lengths {h,h,h) = 

We note that an n-gear has n leaf vertices and n polygon vertices, the latter of which have 
degrees 2, 3, or 4. Let G be an n-gear. We parameterize the edges of G snch that its polygon 
is an oriented cycle with respect to the indnced orientation, and corresponding polygon sides 
and teeth are head-to-head or tail-to-tail. Fignre iTCl indicates the 4 possible neighborhoods 
of a polygon vertex. In order to prodnce pairs as in Fignre [21 we introdnce an anxiliary 
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weight w > 0 and consider the vertex conditions described in FignrelH where oriented edges 
and fnnctions on them are denoted by the same symbol. The corresponding operator is 
self-adjoint with respect to the measnre obtained by weighting the Lebesgne measnres on 
the teeth by w. Namely, if / has restrictions and then we consider 


ll/ll 


2 

w 


(/,/) 


^ pi 

V/ p]{x)+wtl{x)dx 

i=i 


If tp has restrictions and then integration by parts gives 


n pi 

(A/,(p)^ = p''{x)Tri{x)+wti{x)Ti{x)dx 

i=i 

= y] / p[{x)T:[{x)+wt[{x)T[{x)dx-^p'{x)Tii{x)+wt[{x)Ti{x) . 

If / and (f obey the vertex conditions described in Fignre [H then the latter snm vanishes, 
which can be seen by collecting terms vertex-by-vertex and considering the 4 cases shown 
in Fignre [TO separately. In particnlar, {Af,ip)u, = {f,Aip)w, and there exists an {,)w- 
orthonormal basis of eigenfnnctions with eigenvalne seqnence 0 = Aq < Ai < A 2 < .... 

We proceed with the core argnment, which shows that this seqnence is contained in the 
spectrnm of the corresponding operator on the dnal n-gear of G. If the neighborhood of a 
polygon vertex looks like one in the npper half of Fignre [TCl with edges pj and pk, then the 
neighborhood of the corresponding polygon vertex of the dnal n-gear looks like the respective 
one in the lower half, with edges pj and pk, and vice versa. This will allow ns to introdnce 
a linear fnnction, called transplantation, between the spans of non-constant eigenfnnctions 
on G and those on its dnal, which is locally given by 

(3) pi = Pi + wt[ and U = p[ —1[, or vice versa. 


Let / be an eigenfnnction on G with eigenvalne A > 0. Assnme that G has a polygon vertex 
of degree 3 as in the npper right corner of Fignre IlCl Let pj, tj, pk, and tk denote the 
corresponding restrictions of / so that p” = Xpj, t" = Xtj, p'l. = Xpk, t'l. = A4, as well as 

f'(0) = 0, t'kiO) = 0, pj{lj) = = pk{0), and p'j{lj) +wt'j{lj) = p'^{0). 

We show that ([3]) gives rise to a fnnction / that is well-dehned at the vertex pj{lj), obeys the 
desired vertex conditions at pj{lj) and tj{lj), and is an eigenfnnction on pj and tj. Namely, 

(4) p'j = pJ + w t'j' = Xp'j + wX t'j = Xpj. 

Likewise, t” = Xtj, p'l = Xpk, and t'l = Xtk. More interestingly, 

= P-{1,) - tlilj) = X{p,{lj) - t,{lj)) = 0, 

Pji^j) = p'jih) + wt'j{lj) = p 1(0) = :pfc(0) = 4(0), 
p'jilj) = X{pj{lj) +wtj{lj)) = A(1 + '«;)pfc(0) 

= X{pk{0) +wtk{0)) wA(pa:(0) - 4(0)) = pl{0) +wi'f.{0). 

Similar argnments apply to the remaining 3 cases in Fignre ITCl In particnlar, ([3]) gives rise 
to a globally well-defined eigenfnnction / on the dnal n-gear of G. For the sake of simplicity, 
we henceforth assnme that / has restrictions {pi)l^i and (L)jLu and that / has restrictions 
{Pi)i=i and (ti))Li- The general case jnst differs by a redistribntion of tildes. 
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(A) Derived from 3-gears with w = 3/2. 



(B) Derived from 4-gears with w = 2. 


Figure 2. Isospectral quantum graphs with Kirchhoff-Neumann conditions. 


It remains to show that for any eigenvalue A > 0, the transplantation given by ([3]) restricts 
to an injective function between the A-eigenspaces. To this end, we write (E]) in matrix form 


Pj 

ti 


w 


Pi 

i' 


■1 7 V b / ’ V 7 x{i + w)\i -I 

Since every n-gear has A = 0 as a simple eigenvalue, we have shown the following. 

Theorem 1. Mutually dual n-gears are isospectral for every weight w > 0. 

It is worth mentioning that normalizing ([3]) by the factor (A(l -|- leads to a linear 

isometry between the A-eigenspaces. More precisely, we have 

P- = {p\ Pwtff)^ + w{v\ - = (1 + w){p'l + lutf). 

In particular, integration by parts gives 

1 


1 -I- tc 


5 ^/ P>i{^)p>i{^)Pwt!^{x%{x)dx 


2=1 

n 

E 

i=l ' 


XpKx) + W XPi{x)dx + '^p'i{x)pi{x) + wt'i{x)ti{x) = All/ll 

f ^ n 


2=1 


1^- 


As above, the latter sum vanishes since / obeys the desired vertex conditions. Yet, the deh- 
nition ([ 2 ]) has the advantage over its normalized version as it is the same on all eigenspaces. 

Finally, we explain how mutually dual n-gears give rise to pairs as in Figure [2l Roughly 
speaking, each eigenspace decomposes orthogonally under the action of the respective graph’s 
isometry group, and the subspace of invariant elements and its orthogonal complement are 
transplanted separately. We exemplify the method with the help of the graphs in Figure [2A1 
Each of them features a Z 2 -action, given by swapping parallel edges that make up a side of 
its central 3-gon, and a Zs-action, which moves all but these edges. Thus, each eigenspace 
is a unitary representation of Z 2 x Z 3 . Eigenfunctions on which Z 3 acts by multiplication 
by or its square vanish on the central 3-gon. Also, their derivatives on the 3 parallel 

edges of length b sum to zero at the common polygon vertex, equally for length c. For 
such eigenfunctions, the graphs in Figure I2AI essentially reduce to the same 3 subgraphs 
each of which has 2 vertices and 3 edges, and we can tranplant trivially from subgraphs to 
subgraphs. Similarly, eigenfunctions that are odd with respect to the Z 2 -action are supported 
on the central 3-gon, and we can tranplant trivially from 3-gon to 3-gon. Hence, it suffices 
to consider the spaces of Z 2 x Z 3 -invariant A-eigenfunctions on the graphs. However, these 
spaces are linearly isometric to the A-eigenspaces of the 3-gears in Figure [TBI with w = 3/2, 
where the isometry is given by merging the parallel edges of the graphs in Figure [ 2 A 1 
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3. Combinatorial eigenderivative transplantation 

Long before quantum graphs were introduced as model systems in quantum chaos |KS97l 
IKS99b] . they had been studied under different names in chemistry and biology, see |Kuc02] 
and references therein. Notably, |vB85] reduces the spectral analysis of quantum graphs with 
Kirchhoff-Neumann conditions and commensurable edge lengths to that of combinatorial 
graphs with associated row-stochastic matrices. In contrast to |KS99bl rSS06| iBSOTl IOB12] . 
|vB85| rigorously treats the so-called Dirichlet eigenvalues |BK13] . for which ([1]) fails. 

In the context of n-gears, this leads to a combinatorial eigenderivative transplantation. 
We consider pairs of mutually dual n-gears with commensurable edge lengths ■ ■ ■ Jm)- 

Since scaling all edges by some factor leads to eigenvalues scaled by that same factor, we 
may assume that {h,l 2 , ■ ■ ■ ,^m) £ ^+- As is well-known |BK13] . adding or removing a 
vertex of degree 2 carrying Kirchhoff-Neumann conditions leaves the set of eigenfunctions, 
and therefore the spectrum, unchanged. We thus replace each edge of length li by a path 
consisting of k edges of length 1, that is, we subdivide all edges into edges of unit length. 

In the style of |vB85| . we hrst consider Dirichlet eigenvalues A = with fc G Z. 

By virtue of ([2]), the Neumann conditions at the leaf vertices propagate towards the polygons 
so that any A-eigenfunction has vanishing hrst derivatives not only at the leaf vertices but 
also at the polygon vertices in the direction of the respective leaf vertex. In particular, any 
A-eigenfunction is uniquely determined by its restriction to the polygon, which is, in fact, a 
A-eigenfunction of the circle with circumference li + I 2 + ... + Im- On the other hand, any 
A-eigenfunction of this circle gives rise to unique A-eigenfunctions on the n-gears. 

We turn to eigenvalues A for which ([T]) holds. In particular, any A-eigenfunction / on one 
of the subdivided n-gears is determined by its values at the vertices. We differentiate ([T]) 
and set /* = 1 to obtain 

(5) sin(\/A)/'(0) = \/A (^- cos(\/A)/i(0) + /*(!)) . 

If n is a vertex with neighbors P on the polygon and neighbors T on teeth, then the vertex 
condition at v described in Figure [His given by 

(6) cos(\/A)/(n) = ^ I /K) + ^ ^ /K) 

I '\v'&P v'&T 

Thus, the values of / at the vertices give rise to a cos(\/A)-eigenvector of the row-stochastic 
matrix M that corresponds to the random walk on the vertices where edges belonging to 
teeth are taken w times as likely as edges belonging to the polygon. In fact, the A-eigenspace 
is isomorphic to the cos(\/A)-eigenspace of M. For the sake of brevity, we call M the Markov 
matrix of the n-gear, and denote the right-hand side of ([6]) by M[/](n). We note that M is 
irreducible and has period 2 or 1, depending on whether the subdivided n-gear is bipartite 
or not. Moreover, M = D^^AD 2 where Di and D 2 are invertible diagonal matrices, and A 
is the adjacency matrix of this graph. In particular, M is similar to a symmetric matrix and 
therefore has spectrum in [—1,1], where 1 is the simple Perron-Frobenius eigenvalue, and —1 
is an eigenvalue precisely if M has period 2, in which case it is also simple. In the following, 
we give an alternative proof of Theorem [T] Since the eigenvalues of a quantum graph depend 
continuously on its edge lengths |BK13l Theorem 3.1.2], it suffices to show the following. 

Theorem 2. The Markov matrices of mutually dual n-gears with integral edge lengths are 
isospectral for every weight w > 0. 
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Let / be a function on the vertices of a subdivided n-gear. In view of (jS]), we define the 
outward and inward derivatives of / at the vertex v along the edge shared with vertex v' as 

= -M[f]{v) + f{v') and /[„/,„](t) = -/[^y](n) = M[/](n) - f{v'). 

This definition makes any function satisfy the combinatorial version of the desired vertex 
conditions. Namely, if v has neighbors P and T on the polygon and teeth, respectively, then 

v'eP ti'eT Vti'GP ti'eT / 


In particular, if v has the sole neighbor v', then = 0, and if v has degree 2 and 

neighbors v' and v", then = /[!„ For the sake of simplicity, we assume that 

the underlying n-gear has oriented edges and (ti)r=i Figure [TCl the general case 

is obtained by redistributing tildes. Similarly to Section [2l we denote the restrictions of / 
to the corresponding previously-introduced paths by the same symbol, where each of the 
original n polygon vertices appears in as many paths as its degree. We orient these paths as 
their underlying edges, and define the derivative along pi = [tq, Ti, ..., U/J as 


P'iiVj) = 


hj>0. 


similarly for the restrictions This allows to transplant an arbitrary function / via ([3]) 

to obtain a function / with restrictions and on the subdivided dual n-gear, 

which is well-defined by reason of ([7]) and virtually the same argument that showed continuity 
in the quantum graph setting. In order to show that this transplantation maps eigenfunctions 
of M to eigenfunctions of its counterpart M, assume that M[/] = pf for some p G [—1,1]. 
If [n',n,n"] is part of one of the paths (pi)”=i or say p,, then 


P'M + p'iiv”) = -M[f]{v') + f{v) - f{v) + M[f]{v”) 

= p{-f{v') + M[f]{v) - M[f]{v) + f{v")) = 2/ip'(n). 


This entails the combinatorial version of (jl]), namely, M[f]{v) = fif{v) at all vertices v that 
were introduced with the subdivision. We therefore turn to vertices of the underlying n-gear. 
Let V be the leaf vertex with neighbor on the path tj in the lower right corner of Figure [TCl 
and V be the polygon vertex with neighbors Vp and Vt on pj and tj in the upper one. Then 

M[f]{v) = ijivi) = p'jivp) - t'j{vt) = -M[f]{vp) + f{v) - f{v) + M[f]{vt) 

= p{-f{vp) + M[f]{v) - M[f]{v) + f{vt)) = p{p'j{v) - t'j{v)) = fif{v). 

Similarly, let v and v be the polygon vertices with neighbors Vt^, Vp^, vtj,Vpj, and Vp^, on 
4, pk, Pj, tj, pj, and pk, respectively, and Vt^ be the neighbor of the leaf vertex on tk- Then 

(2 + w)M[f]{v) = Pjivpj) + Pk{vp,) +wik{vt,) 

= P'j (% )+Wt'j (Ut,) + Pfc (Up, )+Wtk (Ui, )+wpk (up, )-wt'k {vt ^) 

= Pi-f{Vpj) - wf{vt.) + (1 + w)f{VpJ) + f{v) + wf{v) - (1 + w)f{v) 

= KM[f]{v) - f{vp.) + w{M[f]{v) - fivt^)) + (1 + w){f{vp^) - M[f]{v))) 
= Kp'jiv) +wt'j{v) + (1 + w)pi{v)) 
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Figure 3. Zeta-equivalent simple digraphs. 

The remaining 3 cases in Fignre [TCl follow similarly. Hence, M[f] = /i/ whenever M[f] = p/. 
In order to determine the kernel of the transplantation, we assnme that / = 0. Since tc 7 ^ — 1, 
we have Pi = = 0 for all i. Thus, if [n, v'] is part of one of the paths or (tj)^=i, then 

0 = = -M[f]{v) + /(F) = -pf{v) + /(F), which leads to /(F) = pf{v). 

Hence, f{v) = p!' f{v) for each vertex v on the central polygon of length I = I 1 + I 2 +... + ln- 
If f{v) = 0 for one such vertex n, then / vanishes on the entire polygon, and through 
M[/] = pf on the entire subdivided n-gear. On the other hand, if f{v) F 0 for a vertex v 
on the polygon, then \p\ = 1, that is, p = ±1. Hence, the transplantation is injective on the 
span of yU-eigenvectors with p F ±1- Since subdivided mutually dual n-gears are either both 
bipartite or both non-bipartite, we have proven Theorem [21 In order to derive an explicit 
conjugator for M and M, we note that their ±l-eigenspaces are given by functions that 
satisfy f{v) = ±/(n') whenever v and v' are neighbors. Thus, if M[f] = ±/, then 

f[v,v']i.'^) = -f[v',v]^'^) = -M[f]{v) + /(F) = Tf{v) + /(F) = 0. 

Hence, the ±l-eigenspaces of M are annihilated by the transplantation. Let J_|_ be a rank-1 
matrix that maps the 1-eigenspace of M onto that of M, and annihilates all other eigenspaces. 
If —1 is an eigenvalue of M, we choose J_ analogously, otherwise we dehne J_ = 0. Writing 
M = D^^AD2 as above, the eigenspaces of M are orthogonal with respect to the inner 
product given by the diagonal matrix D 1 D 2 . In particular, we can choose J+ = J 2 nDiD 2 
where J 2 n denotes the 2nx2n all-ones matrix. If T denotes the matrix that corresponds to the 
transplantation, then C = T + J^ + J_ is invertible and satishes MC = CM by construction. 
In addition, the arguments at the end of Section [2] equally apply to the subdivided versions 
of the graphs in Figure [21 which shows that their Markov matrices are isospectral. 

Lastly, we mention further presences of the eigenderivative transplantation method in 
terms of conjugacy. In fact, the method became apparent to us when we discovered cor¬ 
responding conjugators for the matrices H^(\/A) and H^(\/A) in |OB12] . which arise when 
one assumes ([T]) on all edges, and which characterize eigenvalues through the transcendental 
equations det(HF\/A)) = 0. Another characterization of the eigenvalues of a quantum graph 
is given by the scattering approach, which yields an exact trace formula |KS99bl IBE09] . 
It can be shown that mutually dual n-gears have conjugated edge F-matrices, giving yet 
another isospectrality proof for the pairs in Figure [21 

Finally, we consider the digraphs G and G in Figure [3l Note that the teeth of G are 
head-to-tail with their corresponding polygon side. We let I 12 and J 12 denote the 12 x 12 
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identity and all-ones matrix, respectively. For G, we denote its adjacency matrix by Ac-, its 
out-degree matrix by and its in-degree matrix by Dq, the latter two of which have the 
row sums of Aq and on their diagonals. For z = (t, y, a, 13, 7 , 6) G C®, we let 

Lg{^) = xlu + yJi 2 + T Ag + {3 Al^ + 7 + S Dq, 


and similarly for Lq{z). The homogeneous polynomial det(LG( 2 ))|j/=o can be viewed as a 
generalized characteristic polynomial. It determines, and is determined by, the reversing zeta 
function |Herl4] . which generalizes the Bartholdi zeta function |Bar99] to digraphs, which 
in turn generalizes the famous Ihara-Selberg zeta function |Iha 66 ] . The matrix 


T = 


/ 

2 a ^7 

0 

0 

0 

0 

0^(3 

0 

0 

0 

0 

V 0 


a 

Q!^7 

0 

0 

0 

0 

af3‘^ 


0 

0 

2 a ^7 

0 

0 

0 

0 


0 

0 

0 

a^7 


tt/Sy 0^(3 

0 0 

0 0 

0 0 


0 

0 

0 

0 


0 

0 

0 

0 

a^7 


2 a ^7 

0 

0 

0 

0 
3 


a 

0 
0 
0 
0 

/9^7 0 

0 (T/dy a^/d 


0 

0 

0 

0 


0 

0 

0 

0 

0 

—0^(3 

0 

0 

0 

0 

0 


a 

a'^'j 

0 

0 

0 

0 

-af3‘^ 

-af3'y 

0 

0 

0 

3^6 


0 

0 

0 

0 

0 

0 

0 

-0^(3 

0 

0 

0 


0 

0 

0 

a'^'j 

0 

0 

0 

0 

-(3^j 

0 


0 

0 

0 

0 

a^y 

0 

0 

0 

0 

—a/3^ 

-af3'y 


\ 


0 
0 
0 
0 
0 

0 
0 
0 
0 
0 

-a^/d / 


satishes Lq[z)T = TLg{z) and has determinant ((2 q!^)° — (2a^y)“)a®/d^°, which can be seen 
by adding its last 6 columns to its hrst 6 ones. In particular, det(LG(2)) = det(L( 5 ( 2 )), 
meaning G and G are zeta-equivalent and have zeta-equivalent complements [Her 14] . The 
pattern of non-zero entries in T is reminiscent of combinatorial eigenderivative transplan¬ 
tation. The conjugator T can be readily generalized to all pairs coming from dual n-gears 
all of whose polygon vertices have degree 3, where teeth have to be oriented as in Figured 
However, the graphs in Figure [TBI do not lead to zeta-equivalent non-isomorphic digraphs. 
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